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Abstract 

Competition of spatial and temporal instabilities under time delay near the codimension-two 
Turing-Hopf bifurcations is studied in a reaction-diffusion equation. The time delay changes re- 
markably the oscillation frequency, the intrinsic wave vector, and the intensities of both Turing 
and Hopf modes. The application of appropriate time delay can control the competition between 
the Turing and Hopf modes. Analysis shows that individual or both feedbacks can realize the 
control of the transformation between the Turing and Hopf patterns. Two dimensional numerical 
simulations validate the analytical results. 
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I. INTRODUCTION 



Spatiotemporal patterns have been of great interest in a variety of biological, chemical, 
and physical contexts l|. From the viewpoint of dynamics of pattern formation, these pat- 
terns can be classified as three types according to the instability mechanisms: i) periodic 
in space but stationary in time, ii) oscillatory in time but uniform in space (in some cases, 
which often gives rise to phase waves due to the couple of the phases of the oscillators), iii) 
periodic in space and oscillatory in time. Based on linear stability analysis of a homogeneous 
state, these instabilities correspond to Turing, Hopf, and wave bifurcation, respectively. Re- 
cently, the interactions between these instabilities, which often give rise to novel phenomena 



2- 



Yang et al studied the interac- 



of pattern formation, have attracted much attention 
tions between the Turing and wave instabilities and observed the coexistence of Turing spots 
and antispirals {2]. Interaction of two Turing modes leads to the formation of "black-eyes" 
pattern and the superposition patterns combining stripes and/or spots j^. Of particular in- 



;erests are the interaction between the periodic Turing mode and the oscillatory 
4]. It has been explored in a Brusselator modeljs-?!, in a semiconductor model 



iopf mode 



-lioj, and 



in a Lengyel- Epstein model ll|, [12]. The competition of the two modes near codimension- 



two Turing-Hopf bifurcations gives rise to the oscillatory Turing patterns, the stable squares, 
and the bistability between homogeneous oscillatory states and hexagonal Turing patterns. 

Time-delayed-feedback is a practical method for controlling the dynamics of a system 
[isl . \i\ . It is based on applying feedback perturbation proportional to the deviation of the 
current state from the delayed state of the system. It was originally presented by Ott et 



Iq . Recently, time-delayed- 



al to control the chaotic behavior of a deterministic system 
feedback has been used to control the pattern formations in spatially extended systems. 
The time delay with different feedback forms affects the behaviors of pattern formations 



16N24|. It can stabilize the rigid rotation of a spiral or completely destroy a spiral 



much 

16 1 . The tip trajectories of spirals can be controlled by applying time-delayed-feedback 
171 . Il8| . The spontaneous suppression of spiral turbulence by using feedback has been 
studied experimentally in a light-sensitive Belousov-Zhabotinsky reaction and numerically in 
a modified FitzHugh-Nagumo model 19|. Based on a Brusselator reaction-diffusion model. 



Li et al studied the effect of time delay on various pattern formations, such as the traveling 
and standing patterns, the inward and outward spiral waves 2Ch22|. 
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Just recently, Sen et al studied the transition between the Turing and Hopf instabihties 
in a reaction-diffusion model, in which the time delay is incorporated in the kinetic terms 



25|. However, from the viewpoint of experimental realization, it is not always applicable to 
introduce the time-delayed-feedback into the reaction terms directly. In most of the exper- 
imental cases, the time delay acts as an additional feedback and is expressed as additional 
terms in the reaction-diffusion models |131]. This can be found for example in semiconduc- 



tor devices 



26 



27 L nonlinear optical system consisting of a sodium vapor ce 



feedback mirror |28 
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1 and a single 



30|. So, how 



29| . and light-sensitive Belousov-Zhabotinsky reaction 
about the competition between Turing and Hopf instabilities under time delay in a realistic 
system? In this paper, we study the role of time delay played on controlling the competi- 
tion between the Turing and Hopf modes in a model which is originally derived for charge 
transport in a layered semiconductor system ■ We concentrate on the determination 

of feedback intensities on the oscillation frequency, the intrinsic wave vector, and the inten- 
sities of the modes. Results show that time delay with different forms of feedbacks changes 
the dynamical characteristics near the co dimension-two Turing-Hopf bifurcations. Either 
individual or both feedbacks can realize the transformation between the Turing and Hopf 
patterns. The corresponding numerical simulations are also done in two dimensions. 

II. MODEL AND ANALYSIS 

In this work, we study the following two-component react ion- diffusion model: 

11 — 7/ 

^* = 7 ym - + ^"^'^ + ^' (1) 

[v — uY + 1 

vt = aUo -{v- u)) + D,V\ + G, (2) 



here the time delay is applied with the forms: 

F = gMt-r)-u{t)), (3) 
G = gMt-r)-v{t)), (4) 

This dimensionless model is derived from the charge transport in a layered semiconductor 



system SjJ- The variable v represents the voltage across the heterostructure and here acts 
as an activator. The variable u indicates an internal degree of freedom, e.g., the interface 
charge density, and acts as an inhibitor. The parameters Z)„ and describe the diffusion of 



interfaces charge carriers and the induced inhomogeneity of the voltage drop at the interface, 
respectively. T is an internal system parameter describing the tunneling rate, and jo is an 
external current density. The parameter a denotes the ratio of the time scales of u and v. 
The parameters Qu and are the feedback intensities of variables u and f , respectively, r 
is the delayed time. In the analysis, wc first consider a short delayed time, i.e. r is a small 
value. So, we can expand the feedback terms Eqs. (3) and (4) as: 

u(t-r) = u{t)-r^, (5) 

v{t-T) = v{t)-T^. (6) 

So, we obtain: 

(1 + Tgu)ut = 7 r^— - -Tu+ D^V\ (7) 

[V — uy + 1 

(1 + rg^jvt = aOo - {v - u)) + D^V'^v. (8) 

It can be seen that the terms of time is rescaled and determined by the delayed time and 
the feedback intensities. In order to study the competition between the Turing and Hopf 
instabilities under time delay, we first perform the linear stability analysis of Eqs. (7) and 
(8). The uniform steady states of Eqs. (7) and (8) are determined by the intersection points 
of the two nuUclines: 

~ - T« = 0, (9) 



{V - m)2 + 1 

jo - (^; - «) = 0. (10) 

Obviously, the cubic nullcline of Eq. (9) have a form of "N" as indicated in Fig.l. If the 
two nuUclines intersect at a single point lying on the middle branch of the cubic nullcline, 
the system exhibits Turing/Hopf instabilities. If they intersect at one point lying on the 
outer branches of the cubic nullcline, the system owns excitable property. Otherwise they 
intersect at three points, which corresponds to a bistable system. Here, we only concentrate 
on the case of Turing/Hopf instabilities. The uniform steady state of the system reads 

^^0 = jo + Mo- (12) 
Then, we perturb the uniform steady state with small spatiotemporal perturbation {5u, 5v) 




FIG. 1: Nullclines of the model system. -D„=0.1, i:'„=0.348, a=0.08, jo=2.0, T=0.03. 



exp{Xt + ikr), and obtain the following matrix equation for eigenvalues: 

= 0, 




where, 7= 



where, 



• ■ So characteristic equation reads as 



AX'^-BX + C = 0, 

A= {1 + Tg^){l + Tg,), 
B = {1 + Tg^){-a - D,e) + (1 + r(7.)(7 - T - 0^^), 
C={j-T- Duk^){-a - D,e) + a-f. 
The dispersion relations X{k) of the system are defined by the eigenvalues of Eq. (13), 



Ai,: 



2A ■ 



(13) 

(14) 
(15) 
(16) 

(17) 



If the real parts of eigenvalues are positive, the uniform steady state becomes unstable 
to small perturbations. Here, we focus on the temporal instability (Hopf mode) and spatial 
instability (Turing mode), which correspond to k=0 and ky^O, respectively. Figure 2 repre- 
sents the dispersion relation and shows the effect of time delay on the competition between 
Turing and Hopf modes. It can be seen that the application of the time delay changes the 
oscillation frequency of the Hopf instability and the intensities of both the Turing mode 



and the Hopf mode. Therefore, it affects the results of competition between the Hopf and 
Turing modes. In the following we will discuss the effect of time delay on the Hopf and 
Turing instabilities, respectively. 




FIG. 2: Dispersion relation A(A;) with and without time delay. The thick and thin solid lines 
correspond to the real part of the eigenvalue with and without time delay, respectively. The 
thick and thin dash lines correspond to the image part of the eigenvalue with and without time, 
respectively. The feedback parameters are r=0.1, ^^=0.0, §^=5.0. Other parameters are: -D„=0.1, 
L>^=0.348, Q=0.08, jo=2.0, r=0.03. 



For Hopf instability, the most unstable mode occurs at k=0. So the threshold condition 
of Hopf instability is expressed as 

1H = T+ —— . (18) 

1 + rg^ 

and the corresponding critical frequency of temporal oscillation is 



From Eq. (18) it can be seen that the time delay changes the critical value of Hopf 
bifurcation when the feedback and g^ are applied with different intensities. This means 
that time delay with different feedback intensities can introduce the Hopf instability to 
system, which is similar to the cases in refs. [13, 14]. However, application of the feedback 
with identical intensities does not alter the point of Hopf bifurcation, and it only changes 
the oscillatory frequency as indicated by Eq. (19). In Eq. (19), (1 + Tgu){\- + rg^) must be 




FIG. 3: Dependence of the oscillation frequency on the feedback intensities, (a) gu=gv=go, (b) 
Qu^O, gv=0- (c) gu=0, gvT^O. The solid and dash lines correspond to the cases of Re(A)>0 and 
Re(A)<0, respectively. 0^=0.1, D^=0M8, a=0.08, jo=2.0, r=0.03, r=0.1. 



larger than zero in order to make Wq a real value when applying the time-delayed-feedback. 
Figure 3 (a) shows the dependence of oscillatory frequency w on the feedback intensity when 
applying identical intensities of feedback, such that gu=gv=go- It shows that w decreases as 
the feedback intensity increases. Figure 3 (b) and (c) represent the dependence of oscillation 
frequency on the parameter Qu and g^, respectively. The solid lines indicate that the real 
part of Hopf mode is positive. It can be seen that the Hopf bifurcation occurs within a 
range of feedback intensity gu and g^. If an individual feedback g^ is applied, the oscillation 
frequency first increases and then decreases with g^- Within the two critical points 6.2, 
1.2], the Hopf bifurcation occurs. If an individual feedback g^ is applied, the oscillation 
frequency decreases with g^. There are two critical points gy=[3A, 29.9] within which the 
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Hopf bifurcation occurs. Figure 3 (b) and (c) further confirm that the time delay can induce 
Hopf instabihty if gu^dv 

0.6 
0.4 



0.2 
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FIG. 4: Dependence of the intrinsic wave vector on the feedback intensities. Solid line: gu=gv=go, 
dash line: ^^/O, g^=0, dot line: gu=0, gvj^O. 0^=0-1, D^=0M8, a=0.08, jo=2.0, r=0.03, r=0.1. 

For Turing instabihty, we first focus on the critical wave vector /cq and the crit- 
ical value for Turing bifurcation. The bifurcation threshold for Turing instability is 

1 /2 

72-= [r^/^ + (-^j^) Y^"^- At the Turing bifurcation point the critical wave vector is not 

1 /2 

affected by the time delay, and it reads kQ = { ^^ ) . So the time delay does not change 
the critical wave vector and the critical value for Turing bifurcation. This is due to the fact 
that the time delay is equivalent to rescaling time as shown in Eqs. (7) and (8). Beyond the 
Turing bifurcation, the time delay alters the intrinsic wave vector and the intensity of Turing 
mode, i.e. the maximum real part of the eigenvalues as shown in Fig. 2. The dependence 
of the intrinsic wave vector on the feedback intensity is shown in Fig. 4. If an individual 
feedback is applied, such that (/uT^O, gv=0, the wave vector increases monotonously with 
feedback intensity gu as indicated by the dash line in Fig. 4. If an individual feedback g^ 
is applied, such that (7^=0, gvj^O, the wave vector decreases monotonously with feedback 
intensity g^ as indicated by the dot line in Fig. 4. However, if the competition between 
the feedback g^ and gy reaches balance gu=gv=go, the wave vector remains constant while 
changing the feedback intensities, as shown by the solid line in Fig. 4. The three lines 
intersect at one point at the coordinate (^fo, k) = {—1.6, 0.51). 

Near the codimension-two Turing- Hopf bifurcations, the Turing mode and the Hopf mode 
compete. The competetition depends on the intensities of the two modes. Figure 5 shows the 
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FIG. 5: Dependence of the mode intensities on the feedback intensities, (a) gu=gv=9o, (b) Quy^O, 
gv=0- (c) 9^=0, gvT^O. The sohd and dash hnes correspond to the Hopf and Turing modes, 
respectively. L>„=0.1, Dy=0M8, a=0.08, jo=2.0, r=0.03, t=0.1. 



dependence of the most unstable modes (both Turing and Hopf) on the feedback intensities. 
Here, we use the set of parameters in Fig. 1. With this set the Turing mode is stronger than 
the Hopf mode in the absence of time delay. When the feedbacks are applied with identical 
intensities gu=9v=go, the two maximum real part of the eigenvalues which correspond to 
the Hopf and Turing modes decreases with feedback intensities as shown in Fig. 5 (a). The 
Turing mode remains stronger than the Hopf mode. If an individual feedback is applied, 
such that fi^MT^O, gv=0, the intensities of the two modes decrease with gu as shown in Fig. 
5 (b). We want to point out that in this case, the Turing mode is able to become weaker 
than the Hopf mode within the range — 6.0<(/u< — 1.8, which differs greatly from Fig. 5 (a). 
This means that the Turing pattern could transit into Hopf pattern if feedback is applied 
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with appropriate intensity. If an individual feedback is applied, such that Qu—Q, gvj^^, 
the intensities of the two modes increase with as shown in Fig. 5 (c). There exists a 
parameter range 2.2<gy< 15.3, within which the Hopf mode is stronger than the Turing 
mode. 

So it can be seen that by using time-delay-feedback one can realize the transition be- 
tween the spatial Turing pattern and the temporal Hopf pattern, which depends on how the 
intensities of feedbacks are applied. It is practicable to control the self-organized patterns 
by applying individual or both of the feedbacks. As a result of competition between the 
Turing and Hopf modes controlled by time delay as shown in Fig. 5, an original stronger 
Turing mode becomes weaker and then stronger than Hopf mode with increasing individual 
feedback intensity continuously. 

III. TWO-DIMENSIONAL NUMERICAL SIMULATION 

In view of the aforesaid analysis we carry out numerical simulation by using a generalized 
Peaceman-Rachford ADI scheme on a square grid of 100x100 space units. The time step is 
dt=0.02 time unit and the space step is dx—dy— 1.0 space unit. The boundary conditions 
are taken to be no-fiux. 

In order to keep the validation of the Eqs. (7) and (8), we first use a small delayed time 
T=0.1 in Fig. 6. However, extensive numerical simulations have shown that when applying 
feedback with long delayed time and weak intensity, e.g. t=10.0, gu— 0.1, the analysis results 
still remain correct qualitatively because the delayed time and the feedback intensity are 
incorporated into the coefficient of the left-hand terms of the Eqs. (7) and (8). 

Wc concentrate on the competition between the Turing pattern and the Hopf pattern 
by applying time delay. Fig. 6 (a) first shows a coexistence of Turing stripe and Hopf 
oscillation near the codimension-two Turing-Hopf bifurcations in the absent of time delay. 
The time evolution along the marked white dash line in Fig. 6 (a) confirms the stability of 
the coexistence state as shown in Fig. 6 (b). If an individual feedback gu—5.0 is applied, as 
predicted previously in Fig. 5 (b), the Turing mode becomes stronger than the Hopf mode. 
So we obtain pure Turing stripe as shown in Fig. 6 (c). However, if we apply the feedback 
cjf„=— 5.0, the Hopf mode plays a leading role. So we observe oscillatory pattern as shown 
in Fig. 6 (d). 
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FIG. 6: Two-dimensions numerical simulation, (a) coexistence of the Turing stripe and the Hopf 
oscillation in the absence of time delay, (b) time evolution of variable u along the white dash line 
indicated in (a), (c) stripe pattern under time delay, r=0.1, 5„=5.0, ^^=0.0. (d) Hopf oscillation 
under time delay, t=0.1, y„=0.0, 5^=5.0. The other parameters are I?u=0.1, D„=0.348, a=0.08, 
jo=2.0, T=0.03. 

If an individual feedback g^—b.^ is applied, the original Fig. 6 (a) evolves into Hopf 
oscillation as Fig. 6 (d) which shows the dominance of Hopf mode. However, as indicated in 
Fig. 5 (c), if we use gy=~2.5, the Turing mode becomes dominant. We obtain stripe pattern 
like Fig. 6 (c). So the application of the time delay can realize the transformation between 
the Turing and Hopf instabilities resulting from the competition between the Turing and 
Hopf modes. 

IV. CONCLUSIONS 

Wc have studied the competition of the spatial and temporal instabilities under time 
delay near codimension-two Turing-Hopf bifurcations in a two-components reaction-diffusion 
equation. The time-delayed-feedback is determined by the delayed time and the feedback 
intensities. Based on the linear stability analysis, we have investigated the effect of time 
delay on the Turing and Hopf instabilities. Results have shown that appropriate time delay 
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can introduce a Hopf bifurcation when the intensities of feedback are different. It also 
affects the oscillation frequency of Hopf mode. However, the time delay does not change 
the bifurcation threshold and the corresponding critical wave vector for Turing instability. 
Beyond the Turing bifurcation, the time delay begins to work on controlling the intrinsic 
wave vector and the intensity of Turing mode. The wave vector increases (decreases) with the 
feedback intensity when individual feedback gu (or g^) is applied. If the feedback intensities 
gu and g^ are equal, the time delay does not change the intrinsic wave vector. Near the 
codimension-two Turing-Hopf bifurcations the competition between the Turing and Hopf 
modes results in the transition between the Turing and Hopf patterns. The intensities 
of both Turing and Hopf modes decrease (increase) with the intensity of feedback g^ (gv)- 
Analysis shows that individual or both feedbacks can realize the control of the transformation 
between Turing and Hopf patterns. As a result of competition between the Turing and Hopf 
modes controlled by time delay, initially stronger Turing mode becomes weaker and then 
stronger than Hopf mode with increasing individual feedback intensity continuously. We 
have carried out two dimensional numerical simulation which has validated the analytical 
results. Extensive numerical simulations have also shown that when applying feedback 
with long delayed time but weak intensity, the results are still correct qualitatively. We 
hope that our results could be instructive to experiments, such as the charge transport in 
semiconductor devices and the light-sensitive Belousov-Zhabotinsky chemical reactions. 
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